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Abstract 
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by Molecular Dynamics simulations. This information is then introduced into continuum descriptions ac- 
counting for (i) the full curvature and (ii) a long-wave approximation of the curvature (thin film model). A 
comparison of the dependence of the contact angle on droplet size indicates that the theories agree well if 
the contact angles are defined in a compatible manner. 
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I. INTRODUCTION 



In the previous decade increasing attention has focused on the behavior of small quantities of 
liquid on hard [1-7] or soft [8, 9] substrates in equilibrium or under the influence of driving forces 
parallel to the substrate [10-12]. Current research mainly considers two levels of description: 
particle-based models [4, 13-19] and continuum theory [2, 5, 6, 10, 12, 20-22]. The former 
describes the liquid in terms of the position and momenta of particles. These may represent atoms 
in a chemically realistic model or one lumps together a small number of atoms into an effective 
interaction center (called 'bead') in a coarse-grained model. The reduction of the number of 
degrees of freedom and the soft interactions in the coarse-grained description facilitate the study 
of long time and length scales. The properties of particle- based models are studied by discrete 
stochastic simulations, i.e., Monte-Carlo simulation or Molecular Dynamics. The advantage of 
retaining the particle degrees of freedom consists of the ability to refine the model towards a 
chemically realistic description and to include effects of thermal fluctuations and of discreteness 
of matter that are expected to become important on small length scales. However, these stochastic 
simulation techniques are limited to droplets of a linear size that does not exceed a few nanometers. 
Continuum models, in turn, describe the liquid in terms of collective variables that do not refer 
to individual particles. Typical examples of continuum theories are the hydrodynamic description 
in terms of the density and momentum fields or interface models that describe the liquid only 
through a characterization of the motion of its liquid- vapor boundary. Continuum descriptions can 
address engineering time and length scales but depend on phenomenological material constants 
that are often not related in a straightforward way to the microscopic interactions of the particle- 
based description. Thus effort has to be devoted to parameter-passing techniques that transfer 
information from particle-based models to the continuum description. To this end, two questions 
have to be addressed: (i) Which is the relevant information of the particle-based model needed in 
the continuum description and (ii) how can one extract this information from the particle-based 
description in the appropriate continuum form? 

In the present work, we use a coarse-grained particle model of a polymer drop on a solid 
substrate, and a thin film description that characterizes the droplet shape by the location, h, of 
the liquid- vapor interface above the substrate. We explore the behavior of small nano-drops where 
both descriptions are computationally feasible. We extract the interface tensions and the Derjaguin 
or disjoining pressure [23, 24] from Molecular Dynamics simulation of the particle-based model 
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and pass them to contimuum model. Then both approaches are used to determine the equilibrium 
contact angle of a droplet as a function of the size of the droplet and of interaction strength between 
the liquid and the substrate. 

To our knowledge, such a parameter passing scheme has not yet been developped for the case 
of liquid droplets on solid substrates. However, the disjoining pressure itself can be extracted in 
grandcanonical ensemble [16, 25-28]. Additionally, related works exist for other geometries in 
canonical ensemble, such as free standing films or films adsorbed in pores [29] . 

Bhatt et al. [30] extract a disjoining pressure as a function of chemical potential from MD 
simulations for a free standing film of a volatile Lenard- Jones liquid and compare the results with 
the ones of density functional theory. Their approach consists in the definition of the disjoining 
pressure as the difference of normal pressure in the film and the pressure in the homogeneous liquid 
at the same chemical potential as the film. However, as discussed in section III, the measurement 
of the chemical potential in a canonical ensemble is difficult and requires additional simulations. 
Moreover, despite of truncated potentials, they relate the disjoining pressure with solely long- 
range van der Waals dispersion forces and provide therefore comparison to Hamaker theory. The 
short-range forces stay outside the scope of their research. 

A planar liquid film bounded by a solid and vapor is studied by Han [31] using grandcanonical 
MD simulations with a truncated and shifted Lennard- Jones interaction. The disjoining pressure is 
extracted in a similar way as in Ref. [30] and again is associated with only long-range dispersion 
forces. 

Note that parameter passing from MD simulations to continuum hydrodynamics is also fre- 
quently done in the context of liquid flow close to solid substrates [3, 32-36]. However, as 
this works do either not involve free interfaces [34, 36] or do not extract the disjoining pressure 
[3, 33, 35], we do here not discuss them further. 

Our manuscript is structured as follows. In section II we present the particle-based and con- 
tinuum approaches. Then, section III details how we pass the parameters from the particle-based 
model into the continuum description. The subsequent section IV presents the dependence of the 
equilibrium contact angle on droplet size for various interaction energies between the liquid and 
the substrate. In passing, we describe several ways to define the equilibrium contact angle and 
discuss their relation to the macroscopic Young-Laplace law. Section V concludes and gives an 
outlook beyond the case of equilibrium droplets. 
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FIG. 1: Snapshot from MD simulation of a cylindrical drop with illustration of Young's equation 
(left). The enlargement close to the substrate (right) sketches the pairwise bead potentials. 
Coarse-grained beads of polymer chains (blue) interact with each other and with the substrate 
modeled by two layers of face-centered-cubic lattice (lilac) 

H. MODELS 

A. Molecular Dynamics (MD) 

Here, the mesoscopic discrete stochastic description is provided by Molecular Dynamics simu- 
lations of a widely used coarse-grained polymer model [37], i.e., a polymer chain is not represented 
by each and every individual atom but it is modeled as a flexible, linear string of small conglomer- 
ates of atoms. These conglomerates are called "beads". The length of all polymer chains is fixed 
to N p = 10 monomers [11, 38] in all simulations. The potentials used in MD are represented in 
Fig. 1 

All bonded and non-bonded beads have unit mass, m = 1, and interact via truncated and shifted 
Lennard- Jones (LJ) potentials 

U(r) = U u (r) - U u (r c ) (1) 

with 




if their distance is smaller than the cutoff distance r c = 2 x 2 1 / 6 o\ Z7u(r c ) * s me LJ potential 
evaluated at the cutoff distance. All LJ parameters are set to unity, e = 1 and a = 1, i.e., we 
express all energies and lengths in units of e and a, respectively. The reduced time unit r is set by 
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a combination of the LJ parameters as r = a 

The individual beads are connected into chains employing a finite extensible nonlinear elastic 
(FENE) potential given by [39, 40] 

-V 

Ro J 
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FENE 



for r < R Q 

(3) 

oo for r > R 



where R = 1.5a and k = 30e/cr 2 . 

To control the wettability of the polymeric liquid on the substrate we account for the interaction 
of the beads with the solid substrate. The substrate is modeled by a fixed array of atoms as in 
Ref. [11] and not by an ideally smooth and homogeneous wall [15, 38]. Specifically, the substrate 
is represented by two layers of a face-centered-cubic lattice of atoms with a number density of 
p s = 1.5cr~ 3 . We also employ a truncated and shifted LJ interaction between the beads of the 
liquid and the individual constituents of the substrate 

u s (r) = uh(r) - uiAr c ) (4) 

with the length scale cr s = 0.75er. The strength of interaction e s is varied. By changing e s from 
0.2e to e, one tunes the wettability of the system from non- wetting (polymer droplet with a contact 
angle of 6>e = 163°) to complete wetting (polymer film with 9e = 0°). 

All simulations are carried out in a computational domain that corresponds to a three dimen- 
sional box. Periodic boundary conditions are used in the x- and y-directions, whereas the range in 
the ^-direction is limited by a repulsive ideal wall that is positioned far above the polymer liquid. 
The domain side lengths, L x and L y , are chosen in such a way that one may study polymer films, 
L x = L y ), and two-dimensional drops (i.e., ridges in 3d), L y <C L x . These ridges span the simula- 
tion box in y direction and have the cylindrical form whose cross-section is well visible in Fig. 1. 
L y is limited by the Plateau-Ray leigh instability that results in the instability of liquid ridges above 
a critical length. However, as this instability is normally subcritical [41], in a MD simulation L y 
has to be smaller than a critical ridge length L n \ that is smaller than the one resulting from the 
linear stability analysis of a ridge. 

The radius of a 2d drop (3d ridge) scales as y/N (in comparison to N 1 ^ 3 for a spherical 3d 
drop), allowing us to study larger droplets [11]. Moreover, the length of the three-phase contact 
line, 2L y , is independent of the 2d droplet size. Thus, there is no direct effect of the line tension 
on the shape of the droplet. 



The temperature of the system is controlled by a dissipative particle dynamics (DPD) thermo- 
stat [42, 43]. In DPD, the total force on a given monomer is given by 

F tot = ^(F ii + Fg + Fg), (5) 

where the conservative force Fy is derived from the potential between monomer i and monomer 
j, is a dissipative force and F^ is a random force. The dissipative and random forces act on 
pairs of particles and are of the form 

Fg = -7dpd WDO^Xey • v ij )e ij , (6) 

Fg = Cw R (r i ,)%e ij , (7) 

where = |r ; — rj|, and the unit vector e ; j = ry/Vjj points from the j — th to the i— th particle. 
In order to obey the fluctuation-dissipation theorem, the damping coefficient, 7dpd, is connected 
to the amplitude of the noise, (, via the fluctuation-dissipation theorem £ 2 = 2/c B T7 DPD and the 
weight functions are defined as 

| (1 - forr < r c 

<4(r ij )=«> D (r ij ) = \ rJ (8) 

1 for r > r c 

We fix 7 DPD = 0.5 in all our simulations. The term 9y in Eq. (7) is a random noise term such that 
9ij = 9ji and its first and second moments are 

(%} = 0, (9) 
(On (t)e M (f)) = (6^1 + 5 tl 5 jk )5(t - f). (10) 

We use uniformly distributed random numbers [44] with the first and second moments dictated by 
the relations above. 

Since the dissipative and random forces satisfy Newton's third law, they locally conserve mo- 
mentum, i.e., they preserve the hydrodynamics of the flow (in contrast to the dissipative macro- 
scopic behavior in Brownian dynamics). Using this DPD thermostat, we maintain the constant 
temperature, k B T = 1.2e. The equations of motion are integrated with the velocity Verlet algo- 
rithm [45] with a time step At = 0.005t. We performed the simulations on GPU facilities using 
the HOOMD Software [46-48]. 

The MD simulations are used to determine parameters that are passed on to the continuum 
model. Before the parameter passing is described in section III, we introduce in the following 
section the continuum model. 
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B. Continuum model (CM) 



We employ a highly coarse-grained description to characterize the free-energy of a droplet on a 
planar substrate in terms of the position of the solid-liquid and liquid-vapor interfaces. Generally, 
the free energy takes the translationally and rotationally invariant form 

F = 7 sl/ dS + jf dS+ [ dS [ dS'g(\r-r'\) (11) 
Jsl Jly Jim Jsl 

where the integrals extend over the solid-liquid (SL) and liquid- vapor (LV) interfaces [49]. In 
Eq. (11), 7sl is the free energy per unit contact area of the liquid and the solid substrate, and 
7 denotes the liquid-vapor interface tension. The last term of Eq. (11) describes the effective 
interaction between the interfaces, and r and r' are points on the liquid- vapor and solid-liquid in- 
terface, respectively. In the following, we restrict our attention to 2d droplets on a planar substrate 
(cf. Fig. 1), choose the x-coordinate along the planar solid substrate and denote by z = h(x) the 
local distance between a point r = (x, y, z = h(x)) of the liquid- vapor interface and the planar 
substrate (Monge representation). The interaction of a point on the liquid-vapor interface with the 
solid is obtained by integrating over the substrate area 

g(h)= [ dS'g(\r-r'\) (12) 
Jsl 

which for a homogeneous substrate only depends on the distance, h, due to symmetry. g(h) is 
the effective integrated interaction between a point of the liquid- vapor interface with the homoge- 
neous, planar substrate, and it is termed interface potential. In this special case, the free energy 
functional (11) takes the form (up to a constant) 



F[h] = L y J dx^/l + (d x h) 2 [7 + g(h) 



(13) 



where L y denotes the system dimension parallel to the cylinder axis. In the limit that the equilib- 
rium contact angle is small, one can adopt a long-wave approximation (or small-gradient expan- 
sion) 



F[h]^L y J dx [l + ^/i) 2 + ---] [i + g(h) 



(14) 



The equilibrium shape of the droplet is obtained by minimizing this free energy functional subject 
to the constraint of fixed droplet volume 

Kjrop = L y dx h(x) = const (15) 



yielding the condition 



ir(x) 



1 5F 

L y Sh(x) 



A 



(16) 



where A is a Lagrange multiplier constraining the droplet volume. Using Eq. (13) we obtain 



7r(x) 



1 + {d x hf d h g 



9 xx h 



i + g{h) 



d h g 



d x h 



'l + (d x hf 



i + g{h) 



[i + (d x h) 2 f 2 ^i + {d x hf 

In the limit of small contact angles, \d x h\ <C 1, this equation adopts the form 

7r(x) = d xx h 7 + g(h) - d h g. 
The pressure (17) consists of two contributions: the curvature pressure, where «f, 



(17) 



(18) 



ull 



dxx h 



[iMd x hf] 3/2 

is the curvature and 7 + g(h) is the effective tension of the interface a distance h away from the 
solid substrate, respectively. The Derjaguin (or disjoining) pressure H(h) = —dhg(h) models 
wettability [23, 24]. The dimensionless ratio g{h)/^ dictates the shape of a drop in the continuum 
model and it is this parameter that we will extract from the particle-based model in Sec. III. 

A spatially non-uniform pressure, n(x), gives rise to a flow of liquid inside the film. Using the 
Navier-Stokes equation and employing the long-wave approximation [7, 50, 51], one obtains 



d t h = -d x v = -d x {Q{h)d x 7i{x)}. 



(19) 



Here Q(h) = h 3 /3r] is the mobility, r\ is the dynamic viscosity of the liquid. Note, that T is a flux 
that is written as the product of a mobility and a pressure gradient. Eq. (19) with (18) is sometimes 
called a thin-film or lubrication model. 

The equation describing stationary solutions may either be obtained by directly minimizing the 
functional F[h] according to Eq. (16) or, alternatively, one sets d t h = in Eq. (19) and integrates 
twice taking into account that T = in the steady state. Here we use numerical continuation 
techniques [52] to solve the resulting ordinary differential equation as a boundary value problem 
on a domain of size L with boundary conditions such that the center of the resulting drop solution is 
positioned on the right boundary (x = L) and on the left boundary (x = 0) the profile approaches 
a precursor film. The volume is controlled by the integral condition, Eq. 15. Figure 2(a) presents 
typical drop profiles for various volumes whereas Fig. 2(b) gives the maximal drop height as a 
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FIG. 2: (a) Shown are selected half -profiles of droplets at volumes as given in the legend and (b) 
the bifurcation diagram presenting the drop height in dependence of the drop volume. 
Calculations are performed with (i) the full curvature, i.e., Eq. (16) with (17), and (ii) the 
long-wave curvature, i.e., Eq. (16) with (18). Case I and II refer to usage of only 7 or the full 
7 + g(h) as prefactor of curvature, respectively. The profiles in panel (a) are obtained with case I 
for full curvature. The volume is controlled through appropriately adapting the Lagrange 
multiplier A at fixed domain size L = 4000. The employed disjoining pressure and interface 
tensions are extracted from MD simulations at e s = 0.81e (equivalent to an equilibrium contact 
angle of 9e = 23.57°, for details see below section IV). 



function of drop volume. Note, that there exists a minimal droplet volume V sn given by the saddle- 
node bifurcation in Fig. reffig:typical-int(b). If one decrease the volume below V sn , the droplet 
collapses, i.e., it changes discontinuously into a flat film. The transition is hysteretic (first order) as 
the primary bifurcation at V c is subcritical. The situation is different for freely evaporating droplets 
when the chemical potential is controlled instead of volume. For a more detailed comparison of 
the two cases see Ref. [7]. 



III. PARAMETER PASSING BETWEEN PARTICLE-BASED MODEL AND CONTINUUM DE- 
SCRIPTION 

The particle-based model is defined in terms of pairwise interactions between beads, while the 
information that dictates the behavior of the continuum description is the liquid- vapor tension, 7 
and the interface potential, g(h). The latter quantifies the free-energy cost of locating the liquid- 
vapor interface a distance h away from the solid substrate. Several strategies have been proposed 
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to measure the interface potential in computer simulation of particle-based models: (i) The inter- 
action between the interface and the substrate can be obtained in the grandcanonical ensemble, 
where the chemical potential // controls the fluctuating thickness of the wetting layer of the liquid 
on the substrate. The probability, P(h), of observing a wetting layer of thickness h is related to 
the interface potential via g{h) = —k-oThiP{h)+ const [16, 25, 27, 28], where the choice of the 
constant ensures the boundary condition g(h — > oo) = 0. While being elegant, this computational 
technique is limited to simple models because the grandcanonical ensemble requires the insertion 
and deletion of polymers and concomitant Monte-Carlo moves are only efficient for short poly- 
mers, low densities or in the vicinity of the liquid-vapor critical point, (ii) A negative curvature 
of the interface potential at a thickness h signals the spontaneous instability of a wetting layer. 
From the characteristic length scale of this spinodal dewetting pattern one can deduce information 
about d 2 g{h)/dh 2 [53, 54]. (iii) Here we use the pressure tensor. This is a general technique that 
is not limited to short polymers or low densities. It does not require the implementation of parti- 
cle insertion/deletion Monte-Carlo moves and can be straightforwardly implemented in standard 
Molecular Dynamics program packages. 

A. Virial pressure for a liquid film on a solid substrate 

We study a supported polymer film as illustrated in Fig. 3 in the canonical ensemble. By 
virtue of the low vapor pressure of the polymer liquid, one can neglect evaporation effects. The 
fiat liquid-vapor interface allows us to divide the system into thin parallel slabs (separated by the 
horizontal grey lines in Fig. 3), whose normal vector n is perpendicular to the substrate. All 
relevant quantities can then be averaged over each slab, resulting in fields that depend on the 
^-coordinate only. 

In order to obtain the tension of the liquid- vapor and solid-liquid interfaces, 7 and 7 SL , as well 
as the interface potential, g(h), we consider a virtual change of the geometry of the simulation box 
such that the total volume V remains unaltered. Using the scaling parameter A, we relate the new 
linear dimensions, L' x) L' y) L' z of the simulation box to the original ones via L' x = \f\L Xl L' y = 
y/XLy, L' z = jL z . Only the liquid is subjected to this virtual change of the geometry but not the 
solid support. 

The value A < 1 corresponds to a lateral squeezing of the liquid film on top of a solid substrate 
and a concomitant increase of the film thickness h' = jh, where we have assumed that the liquid 
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FIG. 3: Sketch of the slab geometry used to calculate the liquid- vapor interface tension 7. The 
pressure tensor components p n (z) and p t (z) are calculated in every slab k and then their 
difference is integrated across the interface. 



is incompressible. In the continuum model such a transformation gives rise to the the following 
infinitesimal change of the canonical free energy [55] 



dF(A) 



dA 



A=l 



dL'L' 

[7sl + 7 + 9(h)] 



+ 



dg{h) dh' 



A=l 



dh dA 



L X L, 



A=l 



ms dgr(/l), 

1sl + i + g(h) - -^h 



(20) 
(21) 



where, contrary to the related works in grandcanonical ensemble [30, 31], we use the property 
of a canonical one and keep the number of particles in the liquid constant, i.e. constant volume 
hL x L y = h'L' x L' y of the film 

dh' 



dL' x L' y _ 
L' x L' y h! 



0. 



(22) 



The scaling affects the beads of the polymeric liquid only, i.e., the lateral coordinates x and 
y are scaled by the factor y/X and the normal component z is scaled by 1/A. Since the origin of 
the coordinate system is at the center of the simulation cell, the z coordinates of the first layer of 
the substrate is also scaled by the factor 1/A and the second row is shifted such that the distance 
Az between the two layers of the substrate remains unchanged. The lateral coordinates of the 
substrate particles are not changed as depicted in Fig. 4. 
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FIG. 4: Unsealed simulation box (left). Red circles represent two layers of substrate atoms. The 
origin with x and z axes is in the middle of the box. Scaled simulation box (right). The origin has 
moved by ACM against z axis, but its lateral coordinates, x and y, remain unaltered, y axis is 
not shown for simplicity. 



In order to compute the change of free energy, we consider the canonical partition function 



Z 



n!A| n 



i=l 



Yi exp 



(23) 



i<j 



where n is the number of particles in the system, (3 = ^ and At is the thermal de-Broglie 
wavelength. U denotes the bonded and non-bonded interactions between the polymer beads i and 
j, and II s are the interactions between the polymer beads i and the substrate particles s. 

This separation of potentials allows us to express the partition function, Z(\), of the scaled 
system through the scaling transformation of the original positions 



i<j 

-f>Y. u ' 

si,i 



X(xi - xj), \/X(yi - yj), j(zi - Zj) ) (24) 



S2,l 



Axj Xm . V^A?/j Vs\ i \^ Sl 



y/\xi - x S2 , V\yi - y S2 , -Zi - (-z ai - Az] 



where we explicitly separated the interaction of the polymer beads with the first and second layer 
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of the substrate, z' Sl = \z Sl and z' S2 = z' — Az. Differentiation with respect to A yields 
dZ « ,™ <^n,dU dU , dU 



i<j j j j i<j 



/l dU s dU s _ dU s 
^ \ 2 dx isi Xl dy isi Vl dz isi 



31,1 



' E L ( ;>, • - I' 



) exp 


S2,i 


j exp 





2 dxi S2 dyi S2 dzi S2 ^ 

where Xy = x, — x^. Substituting z isi = z iS2 — Az, we write the change of the free energy in the 
form 

dF t T ldZ fin 

dx = ~ k * T zd\ (26) 

^ (l,dU dU dU 

+ (e & Jg* + |> - H*) ) + (g £ ) A ^ (27) 

= ( ^ ] (jz,ij z ij ~ 2 (fx,ij x ij + fy,ijUij) 

+ (E - ^(/u*< + - E •/--, A --) c 28 ) 

\ S,j S2,l / 

where /a,^- denotes the x-component of the force acting between polymer beads, i and j. (• • • ) 
denote averages in the canonical ensemble. 

The first term of Eq. (28) is the anisotropy of the pressure inside the liquid [56, 57]. Using the 
approach of Irving and Kirkwood [58] we define profiles of the normal and tangential pressure in 
a slab k according to 

Pn(k) = k B T(p(k)) + ^;(E ( Vkfaj)), (29) 

and 



V A 

1<J 



Pt(k) = k B T(p(k)) + ^r(j^ ( } (fx,ij x ij + fvjjVij) Vki^ij)^, (30) 

where p(k) is the number density in a slab k and V s i denotes the volume of the slab. The sum 
Yli% mns over particles i and j if the line connecting them crosses the boundary of slab k (then 
r]k(^ij) is the fraction of that line that is located in slab k) or if both particles are in slab k (then 

Tfotry) = 1). 
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Using this definition of the local pressure and Eq. (21), we finally rewrite Eq. (28) as 

Tflim(/i) = 7sl + 7 + g(h) - ~^j^ h = J dz W( z ) - Pt( z )} 




" \ s,i 

The free energy per unit area of the supported polymer film is given by the anisotropy of the 
pressure in the liquid and contributions due to the direct interaction between the liquid and the 
solid substrate. In the limit that the substrate is laterally homogeneous the terms involving the 
lateral forces between solid and liquid vanish. 

We particularly stress that in the canonical ensemble the difference of the film tension 7m m (^) 
and interfacial tensions 7sl and 7 is not the interface potential g{h) [29-31], but of the form of 
Legendre transform g(h) — h^^-. 



B. Solid-liquid and liquid-vapor interface tensions 

In the absence of a solid substrate, the liquid is separated by a liquid-vapor interface from its 
coexisting vapor phase. In this special case, Eq. (31) simplifies and allows us to measure the 
liquid-vapor interface tension through the anisotropy of the pressure tensor components across the 
interface as [5 9-6 1 ] : 

/"2bot 

7= / dz[ Pn (z)- Pt (z)] (32) 

J -Stop 

We find 7 = 0.512 ± 0.006e/cr 2 which agrees well with previous calculations for similar systems 
[1 1]. Mechanical stability requires that the normal component of the pressure is constant through- 
out the system and equals the coexistence pressure [60]. Since the vapor pressure of a polymer 
melt is vanishingly small, p n (z) ~ 0. We also note, that the anisotropy of the pressure is localized 
around the interface and, therefore, the integration can be restricted to an interval [zb ot , -^top] around 
the interface. At the temperature of k^T je = 1.2 the coexistence density of the liquid inside a 
thick polymer film is p a 3 = 0.786. 

If we consider a liquid in contact with the solid substrate, we can measure the solid-liquid 
interface tension 7sl according to Eq. (31). Like in the case of the liquid- vapor interface, the 
anisotropy of the pressure, as well as the additional contribution due to the interaction between 
the liquid and the solid, are localized in a narrow region near the interface between the polymer 
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FIG. 5: The dependence of solid-liquid interface tension 7 S l on the strength of solid liquid 
interaction e s . The horizontal dashed line represents the value of the liquid- vapor interface 
tension 7. 



liquid and the solid. The solid-liquid interface tension depends on the strength e s of the attractive 
interaction between solid and polymer liquid. The simulation results are presented in Fig. 5. 

If the droplet on a substrate depicted in Fig. 1 is at equilibrium, one may describe the equilib- 
rium of forces acting on its contact line by the macroscopic Young-Laplace equation that relates 
the interface energies and the equilibrium contact angle 8 E [62, 63], 

7sl + 7 cos 8 E - 7sv = 0. (33) 

Here, 7 SV is the solid-vapor interface tension. Since the vapor pressure is vanishingly small for 
our polymer melt, we can neglect the interface tension between the solid substrate and the vapor 
phase, 7sv ~ to a first approximation. Using this approximation, we find that the wetting and 
drying transitions occur at 7sL(e s ) w —7 and 7sL(e s ) « 7, respectively. From the data in Fig. 5 
we locate the wetting transition at e s ~ 0.83e and the contact angle reaches 180° for small values 
of e s < 0.2e. 
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FIG. 6: A part of a system used to determine solid-vapor interface tension 7 S v- The droplet serves 
as a reservoir to the chains adsorbed on the substrate. The yellow dotted line indicates the 
curvature of the liquid- vapor interface. The radius of curvature R x is indicated by the orange 
arrow. 





e s = 0.75e 


e s = 0.80e 


e s = 0.81e 


e s = 0.82e 


7sv, [e/o- 2 ] 





-0.00281 


-0.00475 


-0.00523 (-0.01642) 


7sl, [e/ff 2 ] 


-0.32576 


-0.44737 


-0.47419 


-0.49761 


#E (at 7sv = 0), [degree] 


50.50 


29.14 


22.20 


13.69 


6> E , [degree] 


50.50 


29.77 


23.57 


15.98 (20.03) 



TABLE I: Interface tensions of solid-vapor and solid-liquid interfaces and contact angles with 
(6> E ) and without (6> Eo ) taking the solid-vapor interface tension into account. For e s = 0.82e the 
value 7sv is affected by the finite value of Ap and we provide in parentheses an alternative 
estimate of the contact angle. 

C. Solid-vapor interface tension 

While the approximation 7sv ~ is appropriate for small values of the strength of attractive 
solid-liquid interactions, e s , the quality of this approximation deteriorates in the vicinity of the 
wetting transition. If the wetting transition were of second-order, the amount of liquid adsorbed 
onto the substrate, would continuously diverge as we approach the wetting transition. Even for a 

16 



first-order wetting transition we expect that the adsorbed amount will increase when e s increases 
towards its transition value. In this case the approximation 73 v ~ becomes unreliable and we 
employ a meniscus geometry as shown in Fig. 6 to extract the value of the solid-vapor tension. 
The film thickness is chosen sufficiently large, such that the deviation of the pressure from its 
coexistence value, Ap ~ — r +_r > with R x and R y = 00 denoting the principle radii of curvature 
of the meniscus, has only a small influence on the adsorbed amount of polymer and 7sv- Since 
Ap < 0, the adsorbed amount in the simulations will be smaller than at coexistence, 7sv will be 
too large (i.e., negative 7sv will have an absolute value that is too small), and we will slightly 
underestimate the contact angle, 9 E . This correction to the deviation of the approximation 73 v ~ 
0, however, is insignificant for the used system size for all values of e s but the close vicinity 
(e s = 0.82e) of the wetting transition at e s ~ 0.83e, where we have used an alternative methods as 
described in the following Sec. HID. 

For the calculation of 7sv we used the same procedure as earlier for the solid-liquid interface 
tensions of a film, but the procedure is only applied to the part of the simulation box that is far away 
from the meniscus -forming liquid bridge. The values of 7 S v and 7 SL (for comparison) are presented 
in Table (I). One notices the increase in 7 S y when the wetting transition is approached. However, 
compared to the influence on the solid-liquid interface tension the effect is small. Nevertheless, 
it becomes more important the closer one comes to the wetting transition, and the correction of 
the contact angles is significant when one compares profiles of drops of different sizes with the 
prediction of Eq. (33). 

We compare the shape of drops obtained from the particle-based and continuum description in 
the vicinity of the wetting transition. In the following detailed comparison we employ the values 
e s = 0.75e, e s = 0.80e, and e s = 0.81e for the solid-liquid interaction strength that correspond to 
contact angles 8 E = 50.50°, 8 E = 29.77°, 8 E = 23.57°, respectively. For e s = 0.82e, however, we 
will use the more accurate value, 9 E = 20.03°, extrapolated from the interface potential instead. 

D. Interface potential and Derjaguin pressure 

If we consider a polymer film on top of the solid substrate, Eq. (31) provides information about 
the solid-liquid and liquid- vapor interface tensions, 7sl and 7 as well as the interface potential, 
g(h). For a thick film (cf. Fig. 7), the transitions in polymer density at the two interfaces are well 
separated, and the density at the center of the film approaches the bulk coexistence value. In this 
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FIG. 7: Density profiles of a polymer film at e s = 0.80e. The solid black line represents a thick 
film with a bulk region separating solid-liquid and liquid-vapor interfaces. In the case of a thin 
film (dashed red line) it is no longer possible to distinguish the two interfaces. The dotted 
horizontal line indicates the coexistence number density. 



case, also the contributions to Eq. (31) that stem from the two interfaces can be well separated. 
The anisotropy of the pressure tensor at the solid substrate gives 7sl, and the one at the liquid- 
vapor interface gives 7. Thus, the interface potential vanishes, g(h — > 00) — >■ 0, indicating that the 
liquid-vapor interface will not interact with the substrate if the film is sufficiently thick. 

However, upon decreasing the film thickness, the two interfaces start to interact and the con- 
tributions of the solid-liquid and liquid-vapor interfaces can not be separated anymore. The in- 
teraction between the interfaces is quantified by the interface potential, g(h), or equivalently, by 
the Derjaguin pressure H(h) = From Fig. 7 we observe that for small film thickness both 

interface density profiles are distorted, and the density does not reach its coexistence value at the 
center of the film. The distortion of the density profile far away from the interfaces is characterized 
by the bulk correlation length, £ , which therefore sets the length scale of the interface potential 
[64]. 

Since we have determined 7sl and 7 independently, we are able to extract the interface poten- 
tial, g(h), from the simulation data for thin films. To this end, we have to define the location of the 
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liquid-vapor interface, i.e., the film thickness, h. There are several options: Either (i) one deter- 
mines the position where the density equals a predefined value, typically (p liq + p vap )/2 (crossing 
criterion) or (ii) one defines the film thickness via the adsorbed excess (Gibbs dividing surface), 



Ar ads = L x L y J dz [p(z) - p vap ] = [p hq - p vap ] L x L y h 



(34) 



In this work we adopt the integral criterion (34) to define the film thickness. Neglecting the 
vanishingly small vapor density at coexistence, we obtain 



N 

1 ' mon 



^eff 



(35) 



Pliq^4film 

where iV mon is the number of monomers of the liquid inside the simulation box and A film = L x L y 
is the area of the substrate underneath the film. 

We note that both definitions become problematic for film thicknesses where the curvature of 
the interface potential is negative, < 0. In this regime of film thicknesses a laterally extended, 
homogeneous film becomes unstable with respect to spinodal dewetting [7, 53, 65]). However, 
even in this film thickness region, the films can be linearly or even absolutely stable if the lateral 
extension of the simulation box is sufficiently small. The related critical values depend on film 
thickness [66, see, e.g., Fig. 8 of]. In the simulation, we can still obtain meaningful data for 
the interface potential if we restrict the lateral system size to be smaller than the characteristic 
wavelength of the spontaneous rupture process. 

Additionally, we mention that the liquid-vapor interface in our Molecular Dynamics simula- 
tions exhibits local fluctuation of its height (i.e., capillary waves), and the Gibbs dividing surface 
measures the laterally averaged film thickness. The interaction of the liquid- vapor interface with 
«he substrate imparts a lateral correlation length, (, = 2,fj%, onto these interface flnctna- 
tions. These fluctuations give rise to a weak dependence of the interface potential on the lateral 
system size for L x , L y < i.e., the interface potential is renormalized by interface fluctuations. 
Qualitatively, the effect of fluctuations is to extend the range of the potential, i.e., £ = £ (l + u/2) 
withw = M- [67]. 

The interface potential exhibits a minimum at small film thickness, h min . This film thickness 
characterizes the amount of liquid adsorbed on the substrate in contact with the vapor. As illus- 
trated in Sec. IIIC h m i n = except for the close vicinity of the wetting transition. The free energy 
of such a vanishingly thin polymer film is given by 7 fi i m (/v in ) = 7sl + 7 + g(h m ^ ) = 7sv- Thus, 
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the measurement of the different tensions for a planar polymer film provides the value of g(h min ). 

Alternatively, we can use the measured value g(h min ), in turn, to estimate the solid-vapor ten- 
sion, 7 SV . We have employed this strategy for e s = 0.82e, where the finite curvature of the 
meniscus result in a relevant deviation of the pressure from its coexistence value. Extrapolating 
the simulation data to the thickness h m i n « we obtain 7sv = —0.01642. We will use this more 
accurate value, which is not affected by the curvature of the meniscus and that is compatible with 
the interface potential, in the comparison with the continuum model in Sec. IV. 

Since Eq. (31) only provides the Legendre transformation of the interface potential and we 
require an analytical expression for the continuum model, we make an Ansatz for the functional 
form of g(h). Generally, one can distinguish between short-range and long-range contributions to 
the interface potential [64, 68]. The long-range contribution results from dispersion forces between 
the liquid and the substrate. In our particle -based model, however, we do only consider the short- 
range part as our LJ interaction (2) is cut off at r c . Thus, there is no long-range contribution in 
our model in contrast to previous works, when an effective long-range contribution was taken into 
account despite finite interacion cut off [30, 31]. The short-range contribution to g(h) stems from 
the distortion of the interface profile due to the nearby presence of the solid substrate as illustrated 
in Fig. 7, and it is typically expanded in a series of exponentials [64, 68] 

g B (h) = ae- h 'Z - be' 2 ^ + ce" 3 ^ - de~ Ah ^ + . . . , (36) 

In order to obtain g[h) in practice, we fit its Legendre transform g(h) — h-A by a sum of four 
exponential terms like in Eq. (36), and enforce that the interface potential exhibits a minimum at 
h m ™ ~ (there is no precursor film in our MD model) with a value g(h ia i n ), as obtained by the 
measurement of the interface tensions. The resulting fits for g(h) at e s = 0.75e, 0.80e, 0.81e and 
0.82e are given as solid lines in Figs. 8a-8d. The parameters of the fits are presented in Table (II). 

Using the macroscopic Young-Dupre relation, one observes that value of the minimum of g(ti) 
dictates the contact angle [55] 

g(h min ) = 7(cos0 E - 1) (37) 

Much more information can be extracted from the interface potential: (i) The shape of the interface 
potential controls deviations of the drop shape from a spherical cap in the vicinity of the wetting 
transition, (ii) Within the square-gradient approximation the integral of yg~(h) is related to the 
line tension at the three-phase contact line [69]. For all values of e s investigated in the particle- 
based model, the line tension is expected to be negative, (hi) The observation that g(h) increases 
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FIG. 8: Panels (a), (b), (c) and (d) give the interface potential g(h) at e s = 0.75e, e s = 0.80e, 
e s = 0.81e and e s = 0.82e, respectively. They are obtained by fitting the MD results for the 
tension 7mm(^) — 7 — 7sl of films of various small thicknesses (black symbols with error bars) 
by the expression g(h) + hH(h) (dashed black line) obtained employing the first four terms of the 
short-range part of the interface potential g S r(h). The resulting interface potential g{h) is given as 
solid red line. Note, that the minimal value of g min is always reached at vanishingly small 
thicknesses h w Oct, as there is no precursor film in our MD model. 



above zero at intermediate values of h indicates that the wetting transition is of first-order. 



IV. STATIC CASE - SITTING DROPLETS 



In the following we will compare the shape of droplets obtained from the particle-based model 
and the continuum description. This comparison focuses on droplets with small contact angles 
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Parameter 


e s = 0.75e 


e s = 0.80e 


e s = 0.81e 


e s = 0.82e 


a, [e/a 2 ] 


0.13191 


0.06485 


0.05057 


0.06132 


b, [e/a 2 ] 


1.40871 


0.58700 


0.41256 


0.36875 


c, [e/a 2 ] 


1.67606 


0.70902 


0.50249 


0.42963 


d, [e/a 2 ] 


0.58566 


0.25447 


0.18323 


0.15318 




1.51770 


1.26964 


1.14735 


1.00512 



TABLE II: Parameters of the fitting curves of g{h) for the case where the first four terms of the 
short-range contributions [Eq. (36)] are taken into account. Note, that only three parameters are 
independent since there are two additional constraints: The local minimum criterion at h « 0a 
implies d — (a — 2b + c)/4 and the Young-Dupre relation 37 dictates the value g min by setting 
b = a + c — d — (yf min = 1.5a + 0.5c — 2g min . 




x-coordinate, [o] 

FIG. 9: (color online) Profiles of two-dimensional droplets obtained by cutting cylindrical 
droplets obtained in MD simulations [solid noisy line (black online)] for the case e s = 0.82e for 
two values of h max (4.046a and 12.181a). The corresponding spherical cap fit is given as solid 
smooth line (red online). The MD drops are compared with results of the continuum model 
Eq. (16) with the full curvature [Eq. (17)] and in long-wave approximation [Eq. (18)] that are 
given as dashed (green online) and dotted (blue online) lines, respectively. The inset shows a 
zoom into the three-phase contact line region of the smaller droplet. 

< 50° obtained in the particle-based model for strengths of solid-liquid interaction close to the 
wetting transition (e s = 0.75e to 0.82e). Different numbers of polymer chains are used to create 
cylindrical 2d droplets (3d ridges) of varying volumes and hence heights. Data are sampled with 
a frequency of 4000 MD steps. This time interval between two samples corresponds to the Rouse 

22 



25 
20 

a is 

C 
• i— j 

8 io 

o 
o 

N 

5 



" -100 -80 -60 -40 -20 

x-coordinate, [a] 

FIG. 10: Droplet profiles as obtained from continuum theory with full curvature (heavy solid 
lines) for e s = 0.81e and drop height (a) H = 4 and (b) H = 20. Also shown are the spherical 
caps as obtained from the curvature at the drop maxima (heavy dotted lines), and the tangent lines 
at the point of the steepest slope of the profile (thin solid line), and the tangent line of the 
spherical cap profile at precursor height (thin dotted line). Drop height H is defined as difference 
of height at maximum and precursor height. 

relaxation time for a similar polymer liquid tr = 25.6 ± 5r [11]. For small droplets (up to 
600 chains) the sampling lasted 2 x 10 6 steps, whereas for bigger ones (up to 9600 chains) this 
interval was increased up to 10 7 steps, because large fluctuations of the droplet shape occur. As a 
result, every density profile is obtained by averaging over 500 (small drops) to 2500 (large drops) 
snapshots. To extract the droplet shape and measure the contact angle, we use a set of density 
profiles obtained in 10 independent runs. In total, all large droplets are simulated over 10 8 steps. 

The resulting cylindrical droplet snapshots are cut into slices along the invariant y-direction. An 
average over these slices with respect to its center-of-mass position results in the average number 
density profile in the (x, z) plane. A two-dimensional drop profile is extracted by localizing the 
solid-liquid and liquid-vapor interfaces by the crossing criterion for the density as pj nt = (pu q + 
p vap )/2. Examples of profiles are presented in Fig. 9. The resulting profiles are then compared to 
the ones extracted from the employed continuum models, which are also presented in Fig. 9. 
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One popular characteristics of the drop shape is the contact angle, because it is related to the 
balance of interface tensions at the three-phase contact line of a macroscopic drop. For finite- 
sized drops, however, the contact angle is not uniquely defined: (i) One may define a mesoscopic 
contact angle 9 mcs as the slope at the inflection point of the droplet profile. This is often done in 
thin film models [10, 70], however, the steepest slope obtained in this way may not coincide with 
the (larger) macroscopic contact angle even in the limit of large drop size [71]. This corresponds 
to the distinction of macroscopic and microscopic contact angle in Ref. [72]. Moreover, in the 
particle-based model, the inflection point may be located very close to the three-phase contact line 
where liquid-like layering effects of the particle fluid may occur and affect the drop profile [73]. 
(ii) Alternatively, one may define a spherical cap contact angle by approximating the drop profile 
by a spherical cap profile with a minimal radius of curvature R = —1/k, i.e., using the curvature 
at /i max . The resulting contact angle is # sp h = arccos (1 — h max / R) [74]. In the profiles extracted 
from the particle-based model, we extract 6* sp h by only considering the central part of the drop to 
define the curvature. In this way, the calculation is not perturbed by liquid-like layering effects 
or by the short-range interface potential that distorts the liquid-gas interface close to the three- 
phase contact line. The height of the drop is determined as the difference of the highest point 
of the spherical cap and the position of the solid-liquid interface. # sp h converges to the proper 
macroscopic contact angle in the limit of large drop size, but may misrepresent the shape and 
volume of small droplets. 

For the continuum model, the two angles 9 sp h and # mes are illustrated in Fig. 10 that shows 
two droplet profiles h(x) as obtained from Eqs. (16) with (17), their approximated spherical cap 
profiles and the tangents of h(x) at the point of steepest slope (giving 9 mcs ) and of the spherical 
cap profile at the point where it crosses the precursor height (giving 9 sp h). One clearly notes that 
the two measures differ, and that the difference decreases with increasing droplet size. We will see 
below that the two measures do not converge even for very large drops. In the following we focus 
on the spherical cap contact angle 9 sp h- 

The resulting contact angles for drops of various sizes are presented for different e s as open 
square symbols in Fig. 11. Overall, they agree well with the prediction of Eq. (33) that is given 
as horizontal dashed black line (with the standard deviation indicated as a grey hatched region). 
Corresponding results for the contact angle obtained from the continuum model, employing the 
long-wave approximation for the curvature, Eq. (18), and with the full curvature, Eq. (17), are 
given as well. The results for both are shown as solid (case I) and dashed (case II) lines of different 
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FIG. 11: (color online) Contact angles 9 of droplets of different sizes as a function of droplet 
height. Panels (a), (b), (c) and (d) give results at solid-liquid interaction strengths of e s = 0.75e, 
e s = 0.80e, e s = 0.81e and e s = 0.82e, respectively. Square symbols correspond to the contact 
angle measured in MD using a spherical cap approximation of the droplet profile. Dotted and 
dashed thin horizontal lines correspond to the values 6^ and 9 E obtained from the Young 
equation with and without accounting for the measured solid-vapor interface tension 7 S y, 
respectively. Hashed zones show the standard deviation of 9 E . Panel (d) shows additionally as a 
dot-dashed horizontal line the value of 9 E as extracted from the meniscus geometry. The thick 
solid (case I) and dashed (case II) curves (orange, red and blue online) in panels (a) to (d) give the 
spherical cap contact angle and the mesoscopic steepest slope contact angle # mcs obtained 
from the continuum model with full and long-wave curvature [Eq. (16) with Eqs. (17) or (18)], as 
indicated in the legend. Case I and II refer to usage of only 7 or the full 7 + g(h) as prefactor of 
curvature, respectively [cf. Eqs. (18) and (17)]. 
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colors depending on the angle shown (9 mcs or 9 sp h) and the curvature used. Note that both curvature 
models result in identical results for 9 sph because d x h = at the apex of the drop. This is not the 
case for ^ mes . Case I and II refer to the usage of only 7 or the full 7 + as prefactor of curvature, 
respectively [cf. Eqs. (18) and (17)]. 

The angle 6*,^ obtained in the continuum approach agrees well with the result of the MD 
simulations. This is particularly true for case I (only 7 as prefactor of curvature) where sp h 
converges for large drops to the value obtained with the Young-Laplace equation. The deviations 
of case II from case I are small over the entire thickness range for e s = 0.82e, e s = 0.81e and 
e s = 0.80e, but rather large for e s = 0.75e. Note, that 9 mcs does not agree well with the macroscopic 
angle obtained in the MD simulations. In long-wave approximation it is always at least some 
percent smaller than 9 sp h (more so for small droplets). The angle 9 mcs obtained with the full 
curvature differs less from 9 sph , the difference becomes less than one percent for large drops. For 
both curvature models, 9 mcs always decreases monotonically with decreasing drop size. All these 
statements apply for the respective relation between the various curves in case I equally as in case 
II. The various angles calculated in case I are always slightly below the ones obtained in case II. 

Inspecting Fig. 11, one notes a number of further details that warrant to be highlighted: (i) 
A common feature of the particle-based model for e s > 0.80e, shown in Figs, lib - lid, is the 
overshooting of the values of contact angles at thicknesses h ~ 3 — 7a. This effect can also 
be observed in the spherical cap contact angle obtained from the continuum models. It indicates 
that the product of drop height h max and curvature at the drop apex /t max is not a constant any 
more, instead \h max K m£lx \ first increases with increasing volume (before decreasing again), (ii) 
Another detail one notices is the importance of the solid-vapor interface tension, 7 SV , measured in 
Sec IIIC. At e s = 0.75e it equals zero and at e s = 0.80e the macroscopic contact angles are almost 
the same if one neglects 7 S y or properly accounts for it (cf . the dotted and dashed horizontal lines in 
Fig. 1 lb, respectively). However, the difference between the two approaches becomes increasingly 
important with increasing e s , i.e. decreasing contact angle (Figs. 1 lc and 1 Id). Taking a non-zero 
7sv into account becomes crucial close to the wetting transition. There, for rather small values of 
the contact angle (about 15 — 20°) the difference is of the order of 20 — 40% and accounts for 2 — 6°. 
The difference can lead to an incorrect prediction of the contact angle if one assumes 7 S y = in 
the particle-based model. 

Finally, we note that the error bars of the contact angles 9 measured in MD using a spherical 
cap approximation of the droplet profile (open squares in Figs. 11a to lid) are quite large. They 
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increase with decreasing contact angle even in absolute terms. Several possible explanations exist 
for this behavior: (i) In the vicinity of the wetting transition, there are strong capillary waves on 
the surface of the droplet (particularly close to the three phase contact line) [75]. (ii) The crossing 
criterion we apply to define the profile of the drops (po + Pv)/2 is not a unique choice. There 
are other possibilities to define the local interface position based, e.g., on 10-90% or 20-80% rules 
(cf. [56, 76, 77]). 

Next, we compare the drop profiles as obtained from the particle-based model and the contin- 
uum description. For the case of a rather small contact angle, e s = 0.82e, Fig. 9 gives results 
for a very small droplet of /i max = 4.046cx and a larger one with /i max = 12.181o\ The layering 
effects of the particle-based model are rather independent of droplet size. Obviously, the layer- 
ing of the particle-based model is not captured by the continuum model, however, its predictions 
go smoothly through the steps of the profile and always lay between the lateral end points of the 
steps. At the center of the drop, the spherical-cap fit to the particle-based model and the continuum 
results, obtained with Eq. (16) with the full curvature (Eq. (17)) as well as in long-wave approxi- 
mation (Eq. (18)), nicely agree with each other. As cases I and II can not be distinguished by eye 
alone we have only included case I. 

Differences between long-wave and full curvature and the results of the particle-based model 
are only visible in the contact line region. There, the spherical cap is not a good fit to the particle- 
based model. The two continuum models nearly coincide, implying that the long-wave approxi- 
mation for static droplets is still very good for contact angles around 20°. In the contact line region, 
they seem to represent a better approximation to the particle -based model than the spherical cap. 
One should actually expect this, as the continuum models incorporate the Derjaguin pressure as 
measured in the particle-based model. One may conclude that within its limitations the contin- 
uum model describes the profiles rather well if it incorporates the interface tensions and Derjaguin 
pressure from particle -based model. 

The situation differs for larger contact angles as obtained for e s = 0.80e and shown in Fig. 12(a) 
and (b): (i) The deviation from the spherical-cap approximation is more significant than for the 
smaller contact angle and (ii) the continuum model fails to describe the simulation data for the 
smaller droplet size. The difference between the predictions of the different versions of the con- 
tinuum description is small compared to the deviation between the continuum models and the 
particle-based model. Therefore, the reason of the discrepancy is not rooted in the different ap- 
proximations of the curvature. We note that interface fluctuations in a small droplet are strongly 

27 



suppressed. Therefore, one should rather use the bare interface potential than the one extracted 
from a laterally extended film. Since the bare interface potential has a smaller range than the 
renormalized one [67], we expect the profile of a small droplet to be better approximated by a 
spherical-cap shape than that of a large one, which is indeed consistent with the simulation data. 
Out of the same reason, the predictions of the continuum model are more accurate for the larger 
drop than for the smaller one because it uses the renormalized interface potential as input. This ra- 
tional explains why the predictions of the continuum model systematically deviate from the results 
of the particle-based model for small droplet size. For the large droplet, in contrast, the contin- 
uum model succeeds in describing the deviations from the spherical cap shape, which is larger 
than for small contact angles. The profile of the particle-based model lays right in the middle of 
the predictions of the continuum models. The one that fits best is the case I with full curvature. 
Therefore, we conclude that even for contact angles of about 30° all models agree fairly well with 
the particle-based simulations provided the appropriate interface potential is used. 

Finally, we compare the profiles with a rather large contact angle as obtained for e s = 0.75e 
and shown in Fig. 12(c). For the comparison we use a large droplet with h max = 25.494a. The 
difference between the various versions of the continuum models is clearly seen not only at the 
contact line but over the entire droplet profile. The best agreement with the particle-based model 
is achieved for case I with full curvature; all other versions differ more significantly. Therefore, 
we conclude that for contact angles of about 50° only the model with full curvature agrees well 
with the particle-based model, while the long-wave approximation is not valid anymore. It is not 
advisable to apply at E = 50° where it predicts a contact angle 6 mcs that is 20% lower. 

V. CONCLUSION AND OUTLOOK 

The equilibrium properties of polymer droplets have been studied by Molecular Dynamics 
simulation of a coarse-grained particle-based model and a continuum description in terms of an 
effective interface Hamiltonian. We have devised a simple method to compute the interface poten- 
tial for laterally corrugated substrates, which is based on the anisotropy of the pressure inside the 
film. This general computational strategy can be applied to dense liquids of large macromolecules 
and can be implemented in standard Molecular Dynamics programs. Using the so-determined in- 
terface tensions and the interface potential in the continuum model, we find quantitative agreement 
between both descriptions if (i) the full curvature is used in the continuum model for large contact 
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FIG. 12: (color online) Droplet profiles obtained in MD simulations and with continuum models 
are compared for the case e s = 0.80e, apex heights (a) /i max = 6.594a and (b) /i max = 10.469cr 
and (c) e s = 0.75e and apex height /i max = 25.494cx. The solid curves (light blue online) give the 
liquid-vapor interface as obtained in the MD simulation, while the gray solid curves (red online) 
give the corresponding spherical cap fit. Results of the continuum model (CM) Eq. (16) with the 
full curvature (Eq. (17) - green curves) and in long- wave approximation (Eq. (18) - dark blue 
curves) are shown for cases I and II as solid and dashed curves, respectively. For details see main 
text. 



angles and (ii) the size of the drop is larger than the lateral correlation length, £m, of interface fluc- 
tuations. We also find that for contact angles up to about 30 degree the long- wave approximation 
that is normally used in thin film models describes the droplet shapes even quantitatively quite 
well. 

These results demonstrate that the tensions and the interface potential capture the relevant in- 
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formation that needs to be passed on to a continuum model to describe the equilibrium shape of 
droplets, including the deviations from the spherical cap shape in the vicinity of the three-phase 
contact line. This is an excellent starting point for comparing the dynamics of droplets driven by 
external forces, which we will pursue in the future. 
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